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Abstract 

In this paper, we study vacuum bubble collisions with various potentials including gravi- 
tation, assuming spherical, planar, and hyperbolic symmetry. We use numerical calculations 
from double-null formalism. Spherical symmetry can mimic the formation of a black hole via 
multiple bubble collisions. Planar and especially hyperbolic symmetry describes two bubble 
collisions. We study both cases, when two true vacuum regions have the same field value or 
different field values, by varying tensions. For the latter case, we also test symmetric and 
asymmetric bubble collisions, and see details of causal structures. If the colliding energy is 
sufficient, then the vacuum can be destabilized, and it is also demonstrated. This double- null 
formalism can be a complementary approach in the context of bubble collisions. 
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1 Introduction 



Nucleation and subsequent evolution of vacuum bubbles are very interesting issues in gravity, field 
theory, and cosmology. We know that our universe experiences accelerated expansion and this state 
may be described by a positive cosmological constant. In addition, we know that our universe had 
experienced a stage of inflation. The primordial inflation can be described by the existence of a 
false vacuum in a potential of a field and this necessarily requires a phase transition of the field 
from the false vacuum to the current true vacuum. In this context, the nucleation and evolution of 
vacuum bubbles are quite natural to introduce. In addition, such a study of vacuum bubbles can 
be interesting and important. 

The primordial inflation of the current universe will probably be described by a potential with 
some fields that causes inflation [1 . One of the important points of consensus is that the phase 
transition that finishes the primordial inflation seems to be slow-rolling, not the first order phase 
transition [2 . In other words, it is difficult to explain the ending of inflation by the nucleation and 
percolation of true vacuum bubbles, since the inflating universe quickly expands, so that bubbles in 
general cannot percolate sufficiently [3 . Although we know this pessimistic information, the study 
of bubble nucleation and dynamics is still quite important. There are at least two main reasons 
why. 

1. Cosmic landscape and multiverse: String theory requires moduli stabilization and this 
may be possible by flux compactification [4 . The flux compactification can generate a huge 
number of different vacua [5 , the so-called cosmic landscape [6 . If eternal inflation happens, 
then all vacua can be populated by first order phase transitions (quantum tunneling) and each 
bubble will form a pocket universe. The totality of pocket universes in the eternally inflating 
background is called the multiverse. In the multiverse and cosmic landscape, study of the 
nucleation and subsequent dynamics of vacuum bubbles is illuminated again. 

2. In the context of thermal inflation: After the primordial inflation ended, some unneces- 
sary massive particles can be overproduced [7 . Such massive particles should disappear, and 
this may be possible by introducing thermal inflation. In many phenomenological models of 
thermal inflation, it will be ended by the first order phase transition, not the second order 
phase transition [8 . Therefore, our universe may have some signs of the first order phase 
transition, bubble nucleation, and bubble percolation. 

The bubble nucleating process including gravitation was first studied by Coleman and DeLuccia 
[9]. If we include non-minimal coupling effects that can be motivated from string theory, then 
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possible vacuum transitions can be complicated [10 . The subsequent evolution of the bubble 
including gravitation can be approximated by the thin- wall approximation [11]. However, the thin- 
wall approximation itself is an approximation and hence it may lose dynamics at the transition 
region. The related study beyond the thin-wall approximation was studied by numerical method 

mm- 

In this paper, we focus on bubble collision issues. This issue has a long history. First, without 
considering gravitation, there were some studies using analytic methods or numerical methods 
[M]. One of the interesting points of colliding bubbles is that the colliding bubble may induce a 
vacuum transition [15 or a production of particles [16] using the colliding energy. Second, including 
gravitation, it is highly non-trivial to study the dynamics. We may use the thin-wall approximation 
[17] [H]. In the study of Freivogel, Horowitz, and Shenker [19] (for more advanced review, see [20]), 
they discussed two colliding true vacuum bubbles in the de Sitter background: one is flat and the 
other is anti de Sitter. The analysis in itself is very concise and important, but this may not be 
able to describe the dynamics of fields on the colliding walls. For example, this cannot describe the 
vacuum transition behaviors. Hence, we may need further numerical studies. 

There were some numerical studies of bubble collisions with gravitation. Very recently, Johnson, 
Peiris, and Lehner [21 succeeded in studying bubble collisions with gravitation beyond the thin- wall 
approximation. They assumed hyperbolic symmetry from the Birkhoff-like theorem of colliding 
bubbles and assumed initial data from the Coleman-DeLuccia type solutions. They could solve 
Einstein and field equations numerically and observe and report on symmetric/asymmetric bubble 
collisions and vacuum transitions. 

In this paper, our study can be an alternative approach to the problem using the double- null 
formahsm ESES] (note that, [24] studied 5-dimensional anti de Sitter bubble collisions for different 
purposes). We can list possible differences with the previous studies: 

1. Symmetry: When two bubbles are colliding and there are no other perturbations that 
disturb bubble collisions, then we can assume the hyperbolic symmetry. However, if many 
bubbles collide and the background is very complicated, then the Birkhoff-like theorem can 
be highly non-trivial. For the multi-bubble colliding case [18 , the spherical symmetry can 
be a simphfied model; also, it may be useful to try to extend the symmetry to be not only 
hyperbohc, but also planar. 

2. Beyond the Coleman-DeLuccia conditions: At once we fix a certain symmetry (spher- 
ical, hyperbolic, or planar), we can assign various initial conditions for the field within the 
symmetry. These initial conditions can come from the Coleman-DeLuccia solution, but it is 
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not necessary in genera^. The Coleman-DeLuccia solution assumes 0(4) symmetry. There- 
fore, their solutions are de Sitter, Minkowski, or anti de Sitter. Then the initial tension and 
the initial size of the bubble are already fixed. However, in general, bubbles can be massive 
and hence the tension and size can be free parameters [11 , although the probability of the 
bubble can be reduced. It is reasonable that if a massive bubble is nucleated in the de Sitter 
background, then the causal future of the massive bubble can see the effect of the mass, while 
the causal past may not notice such behavior [I2][T3]. In the double-null coordinates, it is 
very easy to implement this kind of massive thick- wall conditions. 

3. Causal structure: The double- null formalism does not have coordinate singularities in gen- 
eral, since it does not use the coordinate time. In addition, the double-null coordinate pre- 
serves causality. Therefore, it is easy to draw causal structures for complicated cases. 

This paper is organized as follows. In Section [2l we discuss our model for bubble collisions in 
the double-null formalism with spherical, planar, and hyperbolic symmetry. In Section [3l we show 
various phenomena of bubble collisions: bubble percolation, dynamics of bubble walls, and vacuum 
destabilization. Finally, in Section HI we summarize our results. 



2 Model for bubble collisions 

In this paper, we study a gravitational system in Einstein gravity with a scalar field that is governed 
by a potential. Then the action S becomes: 



lOTT Z 



(1) 



where R is the Ricci scalar, (j) is the scalar field, and V^(0) is the potential of the scalar field. 
Einstein equations and scalar field equations are as follows: 

G^, = SttT/, (2) 
= (4) 



-"^For a general spherical bubble, it may be difficult to give an exact relation between a spherical slice and a 
hyperbolic slice. However, we can justify the hyperbolic 'approximation' at least for the vicinity of the colliding walls 
when each walls become sufficiently large. Since the lapse function is approximately ~ 1 — 2M/r + 0(r^), around the 
large wall of a spherical bubble, one can approximate de Sitter, flat, or anti de Sitter. In these limits, we can justify 
hyperbolic slices from spherical slices. The planar or hyperbolic symmetry can be the best trial for general spherical 
bubble collisions beyond the Coleman-DeLuccia conditions: we can try to see qualitative results and we can choose 
a certain symmetry by a practical way. The authors thank to a referee to point this out. 
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Figure 1: Left: General bubble percolation. Gray regions are true vacuum regions and white 
regions are false vacuum regions. Middle: Approximately, this can be described by an out-going 
true vacuum bubble and an in-going false vacuum bubble. Right: When two bubbles are colliding, 
we can magnify around the colliding shells, so that the region can be approximated by left-moving 
and right-moving walls. 

In this paper, we will choose one of the following metric ansatz of double-null coordinate: 

1. Spherical symmetry: dS'^ = dO'^ + sin^ Od^p^ 

ds^ = — Qfg (ix, v)dudv + t1{u^ v)dS'^ (5) 

2. Planar symmetry: dE? = dx^ + dy^ 

ds^ = —a^{u^ v)dudv + r'^{u^ v)dR^ (6) 

3. Hyperbolic symmetry: dH^ = dx^ + sinh^ X<^^^ 

ds^ = —a^{u^ v)dudv + v)dH'^ (7) 

Note that 0<6><7r, 0<(/?< 27r, — oo < x^y < oo, and < x < co. 

Figured] shows general motivations of each symmetries. When vacuum bubbles are colliding and 
percolated, there can be a certain true vacuum bubble that is surrounded by a false vacuum region, 
where the outside of the false vacuum region is already transited to the true vacuum (left). We 
may be able to approximate such a situation by an in-going false vacuum bubble and an out-going 
true vacuum bubble with spherical symmetry (middle). In addition, if we focus on the colHsion of 
two vacuum bubbles (red square in left), the colliding region locally looks like two walls colliding 
(right). Such walls can be approximated by the planar symmetry or the hyperbolic symmetry. 



2.1 Comments on mass functions 

To give the proper initial conditions, first we have to understand the value of metric ansatz in the 
static limit. For spherical and hyperbolic cases, we can impose the Birkhoff theorem and obtain 
the concepts of the local mass functions. We know that the general solutions for V{(p) = Vq will 
look like 

ds' = -N{rfdt' + j^ + r'dnl (8) 



where /t: = +1, 0, or — 1, 



iV^ = .-^-^, (9) 



and d^}-^i = dS^ dVL^ = dR^ or dQ-i = dH^ respectively [19]. 
1. The solution in the spherical symmetry is 



dr'^ 

ds' = -N^infdtl + + r'^dS^ (10) 



where 



Here, is the space-like parameter and is the time-like parameter. We define coordinate 
transformation: 

drs = rs^udu^n^ydv, (12) 
a'? f dv du\ 

= f(-— + — )' 

^ \ 's,n 's,v/ 

and obtain the double- null metric ds'^ = —a'^{u^v)dudv + r1{u^v)dS'^ . Note that we can 
choose the in-going direction r^^u < the out-going direction rg^^ > 0. Thus, we can 
show 

N.irs? = (14) 
Therefore, we can identify that the Misner- Sharp mass function is 

m,M = -[l^^^ j. (15) 

2. The solution in the hyperbolic symmetry will look like 



dr'^ 

ds' = -N^{r^fdtl + ^^-j^ + TldH\ (16) 



Ni = -[l + — + (17) 



where 



rh 3 

from the Birkhoff theorem and N'^ is chosen to be negative. Therefore, now th is the space-hke 
parameter and rh is the time-hke parameter. We define coordinate transformation: 

dr\^ — riiudu-\-rhydv, (18) 
a? f dv du\ 

dt^ = -r -— + — ' (19) 



4 \ '^h,n ^h,^>^ 

and obtain the double-nuh metric ds^ = —a\{u^v)dudv + r'^(u^v)dH'^ . Note that we should 
choose both left-moving and right-moving null directions rh,n > and rh,^ > 0. Thus, we can 
show 

Nl = (20) 
Therefore, we can identify that the mass function in the double- null coordinate by 



mh{u,v) = -y (^1 ^^-^ + — ) • (21) 



Note that usual black hole type solutions in Minkowski vacuum will happen for rrih < limit, 
and hence, in this paper, we are interested in mh < 0. 

2.2 Implementation of double-null formalism 

From now, if there is no confusion, then we will omit the subscripts s, p, or h. 
Define 



Ancj) = S (22) 

and use conventions [23] 

h=^, d=^, f = r,u, 9 = r,,, H^ = 5„, Z = S,,. (23) 
a a 

Einstein tensor components are 

Gnu = --{f,u-2fh), (24) 
r 

= ^ (4r/,„ + + 4/5) , (25) 

Gvv = -^{g,v-2gd), (26) 
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where = +1,0,-1 and a = O^x^x spherical, planar, hyperbolic cases, respectively. Energy- 
momentum tensor components are 

Tt = (28) 

Tt = '^V{S\ (29) 

Tt = (30) 

Therefore, simulation equations are as follows: 

/,„ = 2//i-47rrTl, (32) 
5,, = 2gd-AiTrTt, (33) 
f,v=9,u = -K^-^+4nrTi, (34) 

h. = du = -'-^Tt-^-f. (35) 
In addition, we include the scalar field equation: 

Z^ = Wy = - ^ - 7ra^V'{S). (36) 

r r 

2.3 Initial conditions 

We need initial conditions for all functions {a^h^d^r, W, Z) on the initial u = Ui and v = Vi 

surfaces, where we set Ui = Vi = 0. 

Spherical symmetry We have gauge freedom to choose the initial r function. Although all 
constant u and v lines are null, there remains freedom to choose the distances between these null 
lines. Here, we choose r(0, 0) = ro, f{u,0) = r^o, and g{0,v) = r^o, where ruo < and r^o > 
such that the radial function for an in-going observer decreases and that for an out-going observer 
increases. 

It is convenient to choose Vuo = —1/2 and r^o = 1/2; we choose that the mass function on 
Ui = Vi = vanish, since we can think that the initial point is in the causal past of the bubbles. 
Hence, to specify a false vacuum background, for given r(0,0) = ro and 5(0,0) = 6'f, and if the 
field is at the local minimum, then 

„(„,0)=(,-?lIfiM)-"\ (3.) 

Planar symmetry We first choose r(0, 0) = ro, /(0,0) = 0, 5^(0,0) = 0, a{u,0) = 1, and 
a{0,v) = 1, for simplicity. Therefore, h{u,0) = d{0,v) = 0. 
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Hyperbolic symmetry We have gauge freedom to choose the initial r function. Although all 
constant u and v lines are null, there remains freedom to choose the distances between these null 
lines. Here, we choose r(0,0) = ro, f{u^O) = r^^o, and g{0^v) = r^o, where ruo > and r^o > 
such that the function r increases for both of in-going and out-going observers. 

It is convenient to choose r^o = 1/2 and r^o = 1/2; we choose that the mass function on Ui = 
Vi = vanish. Hence, to specify a pure de Sitter background, for given r(0,0) = ro, 5(0,0) = ^f, 
and if the field is at the local minimum, then 

„,0,0)=(l + ?I^f)li)"'\ ,38) 

To compare previous studies, we can choose the initial mass function nonzero. Then, 

.2 o^„\-l/2 



a(o,o)=ri+^^i^+^y 

V 3 ro J 



(39) 



with a free parameter mo- 



Assignments of initial conditions Now it is possible to assign all of initial conditions along 
initial u = u\ and v = v\ surfaces. 

Initial v = V{ surface: We choose 

Si U < iishell, 

S{U, 0) =\ \S,- Si\ G{U) + St iishell <U< ^Ishell + AiX, (40) 
St iishell + A^l < iX, 

where G{u) is a pasting function which goes from 1 to by a smooth way. We choose G{u) 
by 

Tt{u - ^Xshell) 



G{u) = 1 - sin^ 



(41) 



2Aii 

Then, we know W{u^ 0) = S^u{u^ 0). If />: = ±1, then h{u^ 0) is given from Equation (|32]) . since 
f^u = along the in-going null surface. Then, using h{u,0)^ we obtain a{u^O). If hz = 0^ then 
we obtain f{u^O) from Equation (|32]) . 

We need more information to determine d^g, and Z on the v = surface. We obtain d from 
Equation (|35]) . g from Equation (|3l|) . and Z from Equation 

Initial u = u\ surface: We choose 



5(0,^) 



Si V < VsheW, 

\Si-St\G{v)-St ^shell < ^ < ^shell + A^, (42) 
^shell ^ Av <V, 
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where 

'7l{v - Vshell) 



G{v) = 1 — sin 



2Av 



(43) 



li n = ±1, then we obtain d{{)^v) from Equation (|33]), since g^y{{)^v) =0. If z^: = 0, then we 
obtain ^(0,'^) from Equation (|33|) . By integrating d along we have a{{)^v). 

We need more information for /i, /, and W on the = surface. We obtain h from Equa- 
tion (|35]) , / from Equation , and 1^ from Equation (|36|) . This finishes the assignments of 
the initial conditions. 

We used the second order Runge-Kutta method. The details of the numerical scheme and 
convergence/consistency checks are included in Appendix A and B. 

2.4 Simulation parameters 

Finally, we can specify all simulation parameters. In addition, we need information on potentials 
for detailed purposes. 

We illustrate all simulation parameters. 

1. Initial r and Misner- Sharp mass functions: ro, mo 

2. Thickness of the shells: Aii, 

3. Size of the shells: iigheii, '^sheii 

4. Field values: S't, S'f 

5. Potentials 

Especially, we need further comments on the potentials (Figure [2j). 

1. To study bubble percolation, we use the potential 

„,,„=j|^rifi_^ui_A.^, (44, 



where /3 = 0.1 is a free parameter, Vf is the vacuum energy of the false vacuum, and (/)f is the 
field value of the false vacuum. 



2. To study bubble collision with different true vacua, we use the potential 

1\ f 4) 1 



^2 Jcf>, 



(45) 
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Figure 2: Potentials we used. Vf = 0.0001 and or (j)t are chosen to be 0.1. 



where is the vacuum energy of the false vacuum and C2 is a normalization constant with 
^0 

C2 



(46) 



and is the field difference between the false vacuum and the true vacuum. 

3. If we tilt the potential V2 as follows: 

where Ds/Vf = 0.5 is fixed, we then observe the asymmetric bubble collisions. 

4. To study vacuum destabilization by bubble collisions, we use the potential: 



(47) 



Vf 

04 



d4> + Di [^] , (48) 



where Vf is a constant, D4/V1 = 1.5 is fixed, and C4 is a normalization constant with 



a 



-f 



d(j). 



(49) 



Note that in this case, if we choose a certain S'f , the local maximum is shifted to 0.9502 x S'f, 
and we considered this effect. 



Therefore, for all cases, Vf and Sf = v47r^f (or S't = v47r^t) are free parameters of each potential. 
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ro = 50, Vi, spherical symmetry 
Sf=0.1 Sf=0.2 




Figure 3: Bubble percolation with spherical symmetry: Vi with Vi = 0.00001, tq = 50, Au = Av = 
20, iXsheii = ^sheii = 30, = 0, and varying S'f = 0.1, 0.2, 0.3, 0.4. 

3 Various phenomena of bubble collisions 
3.1 Bubble percolation 

First, we study bubble percolation. By percolation, we mean that two true vacuum bubbles are 
collided and have eventually emerged to one true vacuum region. The question is whether this 
process can continuously happen, and what is the new phenomena that could not be observed by 
the thin-wall approximation? 

3.1.1 Spherical symmetry 

Figure [3] shows the collision of an in-going false vacuum bubble and an out-going true vacuum 
bubble. In this case, we choose the same field value for the true vacuum. Therefore, after the 
collision, bubbles should be emersed and be percolated. In these calculations, we used the potential 
Vi with 0.00001. We also used ro 50, Au = Av = 20, iisheii ^sheii 30, St 0, and 

varying = 0.1,0.2,0.3,0.4. 
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Figure 4: A schematic diagram of bubble percolation with the spherical symmetry. 



As /Sf increases, the tension and the total energy increase. Therefore, the location of the apparent 
horizon and singularity becomes lower and lower as increases; this implies that the apparent 
horizon increases, and a singularity relatively quickly appears. 

In Figure [3l the yellow region is the false vacuum region, while the skyblue region is the true 
vacuum region. The black colored region is around the local maximum of the potential Vi and 
hence it shows the dynamics of the thick wall. 

The behavior of the out-going true vacuum wall before colliding (for v < 30) is changed as 
varies. If is sufficiently small, then initially it shrinks (upper left). However, as Sf increases, it 
tends to expand (lower right). Unless the shell has sufficient energy, it tends to collapse. In addition, 
if the shell has sufficient energy, it can expand, although the inner region becomes destabilized; the 
inside of the out-going shell {v < 30 and u > 50) is not exactly zero and the fields are slightly 
perturbed (we can compare the upper left and upper right diagrams). These results are consistent 
with that of the false vacuum bubble cases [T2\ . 

The upper left of Figure [3] shows interesting behaviors: after the collision, some fields can roll 
up to the false vacuum region (50 < v < 60 and u > 80). This is due to the fact that the colliding 
energy perturbs and pushes the field to the other local minimum. This is related to the vacuum 
destabilization by bubble collisions, and it will be discussed further in Section [331 

We can intuitively draw as in Figure HI Initially, there is an in-going false vacuum bubble 
and an out-going true vacuum bubble (left). The false vacuum region becomes narrower (middle). 
Eventually, the collided region collapses and shrinks beyond the apparent horizon (right). This can 
be compared to the Figure [16] in Section 13.31 This mimics the situation that bubble collisions can 
generate a black hole: Figure [5] 
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singularity 




Figure 5: Causal structure of percolating bubbles with spherical symmetry. 



ro = l^ Vi^ planar symmetry 



Sf=0.1 



Sf=0.2 




20 40 60 80 

V 

Sf=0.3 



20 40 60 80 

V 

Sf=0.4 




20 40 60 80 



20 40 60 80 



Figure 6: Bubble percolation with planar symmetry: Vi with = 0.0001, ro = 1, A^^ = Av = 20, 
'^sheii = '^sheii = 30, = 0, and varying = 0.1,0.2,0.3,0.4. Black curves are contour curves of r 
and the differences of each of the curves are 1 and drawn from r = to r = 100. 
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3.1.2 Planar symmetry 

Figure [6] shows the bubble percolation of two planar bubbles. We assume these conditions: Vi with 
= 0.0001, ro = 1, Au = Av = 20, iisheii = ^sheii = 30, 5't = 0, and varying = 0.1, 0.2, 0.3, 0.4. 
The first clear observation is that after the collision, walls oscillates back-and- forth (upper left). 
As a wall oscillates, oscillating field amplitudes decreases (upper right) and eventually the wall 
disappears (lower left). 

When the tension of the shell is sufficiently small, the shell initially expands to the true vacuum 
region, although it eventually and quickly shrinks to the center (upper left). As the tension of the 
shell increases, the initial expansion of the shell becomes slower and slower (upper right and lower 
left). Eventually, if the shell energy is sufficiently large, it slowly approaches and it may not be able 
to colhde, since the false vacuum region inffates too quickly (lower right). 

We need more comments on the lower right of Figure [6l For all diagrams in Figure [6l we can see 
the details of the r contours. In this case, r is not a real radius; rather, this is just a metric function 
that relates the planar direction and the orthogonal direction. For the inside of the false vacuum 
region, it is space-like, while it is time-like for the true vacuum region. The r function increases 
from lower left to the upper right; and it diverges in the lower right of Figure [6l If we only see 
this diagram, it is not easy to interpret the space-like piece in this diagram, where the numerical 
simulation cannot work. It seems to not be a singularity, since two bubbles are not colHded yet. 
To interpret correctly, we have to compare Figure [71 in Section 13.1.3) in the hyperboHc symmetry, 
we can interpret that the region is future infinity, where the time-like parameter diverges. Then 
we can finally interpret the figure that two bubbles cannot collide since the false vacuum region is 
rapidly expanded. 

3.1.3 Hyperbolic symmetry 

Massless case In Figure [71 we observe bubble percolation with the hyperbolic symmetry. The 
results are similar to the planar symmetric case. We used Vi with = 0.0001, ro = 1, Au = Av = 
20, UsheW = ^sheii = 30, St = 0, and varying = 0.1,0.2,0.3,0.4. 

When two bubbles collide, the collided region oscillates back-and- forth (upper left). The collided 
region eventually decays and disappears (upper right). Similar to the planar symmetric case, as the 
tension increases, the wall slowly moves. For all cases, r contours are always space-like. This is not 
strange, since our metric ansatz assumed r as a time- like parameter. Therefore, r works as a global 
time function. For figures in the lower left and lower right in Figure [71 the function r diverges as it 
approaches the upper right region of the figure. Therefore, we can interpret that two bubbles may 
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ro = l, Vi, hyperbolic symmetry 
Sf=0.1 Sf=0.2 




Figure 7: Bubble percolation with hyperbolic symmetry: Vi with = 0.0001, ro = 1, Au = Av = 
20, UsheW = '^sheii = 30, St = 0, and varying Sf = 0.1,0.2,0.3,0.4. Black curves are contour curves 
of r and the differences of each of the curves are 10 and drawn from r = to r = 1000. 

not be able to collide and the space-like boundary is the future infinity of the false vacuum region, 
where the in-going and out-going null boundaries are future infinity of the Minkowski regions. 

To summarize, we obtain two interesting causal structures in Figure [HI left is normal bubble 
percolation when the tension is sufficiently small; right is the large tension case and due to the 
repulsiveness of large tension walls, bubbles cannot collide and bubble walls touch the future infinity. 

Massive case We can also choose the mo < case. In Figure [H we especially choose — mo = 
ro = 100, since mo and ro should be the same order to see the effect of the mass term. We used Vi 
with Vf = 0.0001, Au = Av = 20, Usheii = ^sheii = 30, = 0, and varying = 0.1, 0.2, 0.3, 0.4. 

As 6'f increases, it is qualitatively similar to the massless case; however, the wall collides relatively 
rapidly than the massless case. Therefore, the existence of mass allows a wide range of tensions 
that allows collision. Therefore, for convenience, we choose this kind of initial conditions for the 
latter sections, since we want to build colliding bubbles. 



17 




planar/hyperbolic symmetry, planar/hyperbolic symmetry, 

bubble percolation bubbles cannot collide 



Figure 8: Causal structure of percolating bubbles with planar/hyperbolic symmetry. Left: bubbles 
are percolated. Right: when the tension increases, bubbles may not be able to collide. 



ro = -mo = 100, Vi, hyperbolic symmetry 
Sf=0.1 Sf=0.2 




Figure 9: Bubble percolation with hyperbolic symmetry: Vi with = 0.0001, tq = —mo = 100, 
Au = Av = 20, UsheW = '^sheii = 30, St = 0, and varying Sf = 0.1,0.2,0.3,0.4. Black curves are 
contour curves of r and the differences of each of the curves are 10 and drawn from r = 100 to 
r = 1100. 
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ro=-mo = 100^ Vi, hyperbolic symmetry 



St=0.1 St=0.2 




Figure 10: Bubble collisions of symmetric different true vacua with hyperbolic symmetry: V2 with 
Vi = 0.0001, ro = —mo = 100, /S.u = /S.v = 20, 'Usheii = '^sheii = 30, 5'f = 0, and varying 
=0.1, 0.2, 0.3, 0.4. Black curves are contour curves of r and the differences of each of the curves 
are 10 and drawn from r = 100 to r = 1100. 




planar/hyperbolic symmetry, 
bubble collision 



Figure 11: Causal structure of colliding symmetric bubbles with different field values. 
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3.2 Dynamics of bubble walls 

Second, we study dynamics of bubble walls. In this case, we especially focus on the case that the 
left true vacuum region and the right true vacuum region have different field values. Therefore, the 
bubble wall region will not disappear. In this section, we will observe their back-reaction to the 
causal structure. Also, we see when the initial conditions of the left and the right are asymmetric. 

3.2.1 Symmetric bubbles 

In Figure [TOl we studied bubble collisions of symmetric bubbles with different field values of true 
vacuum region: V2 with = 0.0001, ro = —mo = 100, Au = Av = 20, UsheW = '^sheii = 30, 6'f = 0, 
and varying = 0.1,0.2,0.3,0.4. Note that we only used hyperbolic symmetry, since the planar 
case is qualitatively similar to the hyperbolic case. 

As St increases, the oscillating frequency of the thick wall increases and is quickly stabilized 
for sufficiently large St. One interesting thing to notice is that after a sufficient time, two regions 
cannot be emersed, since two regions have different field values. Therefore, between two true vacuum 
regions, there should be always a false vacuum region (blue colored region). In a finite advanced 
and retarded time, the function r diverges. Therefore, we can interpret that the coordinate time 
diverges and hence there is the future infinity. For the true vacuum part, the future infinities are 
null. Between two future null infinities, the false vacuum region can roll a crossing point between 
two null boundaries. This is a quite different behavior than the bubble percolation case (Figure [TT]). 

3.2.2 Asymmetric bubbles 

In Figure [121 we studied bubble collisions of asymmetric bubbles with different field values of true 
vacuum region. First we change tensions of two walls: V2 with Vf = 0.0001, ro = —mo = 100, 
Au = 20, Av = 40, UsheW = '^sheii = 30, 5'f = 0, and varying S't = 0.1,0.3. Second we change 
vacuum energy between two true vacuum regions: Vs with Vf = 0.0001, ro = 100, Au = Av = 20, 
^sheii = ^sheii = 30, 6'f = 0, and varying 6't = 0.1, 0.3. 

Asymmetry of tension If one side has larger tension (in our case, the skyblue side has larger 
tension, since the thickness of the shell is narrower than the yellow side), it can push the smaller 
tension region, since it is related to the energy of the shell. The upper two figures show this behavior. 

Asymmetry of vacuum energy If one true vacuum side has smaller vacuum energy (in our 
case, the yellow region is anti de Sitter and hence it is smaller than the skyblue region), it can 
push the larger vacuum energy region, since the difference of the vacuum energy between the false 



20 



ro = -mo = 100^ Vi, hyperbolic symmetry 



St=0.1 St=0.3 




Figure 12: Upper: Bubble collisions of asymmetric different true vacua with hyperbolic symmetry: 
V2 with Vi 0.0001, ro = -mo = 100, /^u = 20, Av 40, u^hew VsheW 30, 6'f = 0, and 
varying S't = 0.1,0.3. Lower: Bubble collisions of asymmetric different true vacua with hyperbolic 
symmetry: V3 with = 0.0001, ro = —mo = 100, Au = Av = 20, UsheW = '^sheii = 30, = 0, and 
varying 5't = 0.1,0.3. 
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i() = -m{)= 1 00, V:;, hyperbolic symmetry 
St.=0.3 




40 80 120 160 



y 

Figure 13: Bubble collisions of asymmetric different true vacua with hyperbolic symmetry: Vs with 
Vf 0.0001, ro -mo 100, Au = Av = 20, ^Xsheii ^sheii = 30, 6'f = 0, and St = 0.3. Black 
contours are r and the differences of each of the contours are 10 and drawn from r = to r = 1000. 
Red and blue curves are v y — and v u — contours. Black arrows are directions where r increases. 

vacuum and the true vacuum is changed to the pressure of the wall. Therefore, in the lower figures, 
the yellow region pushes the skyblue region. 

The lower right in Figure [12] is important to emphasize. This shows the collision between a de 
Sitter vacuum and an anti de Sitter vacuum. This problem was systematically discussed by fi9\ . 
However, here we can see details from numerical calculations. 

Figure [13] is the detailed version of the lower right of Figure [12] Black contours are r and the 
distances of each of the contours are 10. Red and blue curves are r^y = and r^u — contours. 
Black arrows are directions where r increases. Therefore, for the de Sitter region, r contours are 
always space-like. However, inside of the anti de Sitter region, the signature of r changed. After r 
changes the signature two times (r < and r^y < 0), eventually r decreases to zero. This whole 
smooth evolution of the bubble collision and the formation of big crunch singularity are possible 
due to the numerical calculations. 

We summarized in Figure [H] This is the advanced version of the previous work [19 . The left 
bubble is de Sitter (in principle, Minkowski is possible). The left bubble is collided into the right 
anti de Sitter bubble. The wall moves in a time-like direction, is biased to the left direction, and 
eventually becomes parallel to the null direction; the tension via the difference of vacuum energy 
is larger for the anti de Sitter side and the speed of the wall will eventually approach the speed 
of light. The future infinity of the left side is space-like (for Minkowski, it should be null). There 
is a future infinity on the wall, and the wall will touch the future infinity of the left side. When 
the wall touches the future infinity of the left side, the wall may or may not generate a Cauchy 
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Figure 14: Causal structure of bubble collisions between de Sitter and anti de Sitter bubbles. We 
denote local horizons and Bousso wedges. 

horizon; this depends on initial conditions. However, at least in the data of Figure [131 we can see 
that there should be a Cauchy horizon in the anti de Sitter side (green dotted line) because of two 
reasons: (1) the wall almost moves in the left-going null direction and (2) the right side of the wall 
is approximately pure anti de Sitter as time goes on, and therefore, the blue curve {r^u = horizon) 
should be parallel to the null direction. Therefore, the wall and the r^u = horizon will not meet 
each other, and hence, there should be a Cauchy horizon for the anti de Sitter side. We also notice 
the Bousso wedges, that implies the direction where r increases. Due to the horizons, the Bousso 
wedges are changed. 

There are two main differences between the thin-wall approximation [19] and our numerical 
studies: 

1. The crossing point of two horizons {r^y = and r^u = 0) can happen before the causal future 
of the bubble collision, while this was not allowed by the thin-wall approximation. This is not 
strange, because the horizons are dynamically deformed on the wall and connected between 
de Sitter and anti de Sitter regions. 

2. The future infinity can have the Cauchy horizon as we have commented. However, the time 
parameter r will diverge on the wall when the wall touches the future boundary. Therefore, 
in terms of the analytic coordinates of the anti de Sitter side, they should be identified to the 
time-like r = oo boundary (gray diagrams beyond the Cauchy horizon in Figure [H]), although 
we cannot determine beyond the Cauchy horizon from the past data. 
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ro = 50, V4, spherical symmetry ro = -mo = 100, V4, hyperbolic symmetry 

Sf=0.1 Sf=0.1 




Figure 15: Transition to deeper vacuum. Upper: Spherical symmetry: V4 with Vf = 0.00001, 
To = 50, Au = Av = 20, i^sheii = '^sheii = 30, 5't = 0, and 5'f = 0.1. Lower: Bubble collisions of 
asymmetric different true vacua with hyperbolic symmetry: V4 with V[ = 0.0001, ro = —mo — 100, 
Au — Av — 20, i^sheii = '^sheii = 30, 5't = 0, and = 0.1. 
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Figure 16: Deeper vacuum transition by bubble collisions, for spherical symmetric case. 



3.3 Vacuum destabilization to deeper vacuum 

Finally, we focus on the vacuum destabilization to deeper vacuum by bubble collisions. We tested 
two cases. First, for spherical symmetry, V4 with Vi = 0.00001, ro = 50, Au = Av = 20, UsheW = 
^sheii = 30, 5't = 0, and 5'f = 0.1. Second, for hyperbolic symmetry, V4 with Vf = 0.0001, 
ro = —mo = 100, Au = Av = 20, Usheii = '^sheii = 30, = 0, and = 0.1. 

Figure [15] shows the spherical and hyperbolic cases. The skyblue region is the deeper vacuum. 

For the spherical case (left of Figure [T5|). after two walls are collided, a deeper region is generated. 
The boundary between the deeper vacuum and the true vacuum moves along the out-going null 
direction. Hence, it can be outside of the event horizon, if we properly choose the tension. This can 
be interpreted by Figure [161 After the bubble collision, the deeper vacuum region expanded along 
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spherical symmetry, 
vacuum destabilization 




planar/hyperbolic symmetry, 
vacuum destabilization 



Figure 17: Causal structure of deeper vacuum transitions via bubble collisions with spherical (left) 
and hyperbolic (right) symmetry. 

the out-going direction, and hence it can be outside of the apparent horizon. 

For the hyperbolic case (right), the deeper vacuum region and the big crunch singularity of the 
anti de Sitter region are symmetrically generated, as one expected. To summarize, we draw the 
causal structures for vacuum destabilization for spherical and hyperbolic cases in Figure [T71 

4 Discussion 

In this paper, we studied bubble collisions with gravitation by numerical calculations. We used 
the double- null formalism to implement numerical simulations. We used three types of symmetries: 
spherical, planar, and hyperbolic. Three kinds of symmetries can be qualitatively useful in various 
situations during the first order phase transition (Figured]). In this discussion, we will summarize 
our results. 

First, we tested bubble percolation by bubble collisions. After bubbles collided, the false vacuum 
region disappeared and two true vacuum regions are emersed. For the spherical symmetric case, 
it will eventually induce a black hole in a true vacuum background (Figure [5]). For the planar 
or hyperbohc symmetric case, two bubbles are emersed (left of Figure [8]). However, as tension 
increases, if the wall is not sufficient to push the false vacuum region, then two walls may not be 
able to meet (right of Figure [8]) . 

Second, we tested bubble collisions with different field values to see the dynamics of bubble 
walls. For symmetric cases, the collided wall moves symmetrically and the wall will touch the 
future infinity between two true vacuum regions (Figure [TT]). For asymmetric cases, the wall is 
biased due to the difference of pressures. The interesting example is the collision of the de Sitter 
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vacuum and the anti de Sitter vacuum, and the causal structure is analyzed in Figure [Ml 

Finally, we tested vacuum destabilization due to the bubble collisions. For the spherical sym- 
metric case, the deeper vacuum region can be outside of the event horizon (left of Figure [T7|). When 
we compare this to Figure [H we can notice that this can give a significant difference for outside 
observers. For the hyperbolic symmetric case, the deeper vacuum region is expanded over the null 
directions and there will be a space-like big crunch singularity (right of Figure [T7|). 

There were a number of previous studies on bubble collisions, and this paper can confirm these 
studies. In addition, there are mainly new contributions provided by this paper. 

• The spherical symmetric bubble collision (including black hole formation) was studied numer- 
ically with gravitation (Figure [3] and left of Figure [T5|). This can give a general intuition for 
multiple bubble collisions (Figures Hand [my 

• We observed more general initial conditions than Coleman-DeLuccia type bubbles. This was 
easily done, because we used the double-null formalism. One interesting observation is that 
as tension increases walls slowly moves and they may not be colHded (lower right in Figure [6] 
and lower left and right in Figure [7]) • 

• The smooth transitions of metric functions (especially r) are realized. This cannot be done 
by the thin-wall approximation. We obtained future boundary formation from smooth initial 
data without assuming thin- wall and classified various causal structures. Moreover, Figure [T4l 
shows more complicated behavior of horizons of r during bubble collisions. 

There can be further applications and there will be many other models and issues for bubble 
collisions. Also, there can be many problems that can be solved by double- null formalism with 
hyperbolic or planar symmetry. We leave these for future work. 
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Figure 18: Consistency test for Equation (|66l) along = 20, 40, 60, 80. 
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Appendix A: Numerical scheme 

The integration scheme that we used is the second order Runge-Kutta method. Here, we briefly 
discuss the integration scheme. As we see in Section [2^ we organized all equations by a set of flrst 
order differential equations. Let us assume that we know all data for initial surfaces u = Ui and 
V = diS we discussed in Section [231 and observe how we can determine the point {u = Au-\-u^^ v = 
Av + Vi). For example, to obtain the function r, we have to solve 

r,u = /, (50) 
f^u = Right hand side of Equation (|32]) . (51) 

Of course, this is not the unique choice and we can use 

r,v = 9. (52) 
g^v = Right hand side of Equation (|33|) . (53) 
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Figure 19: Convergence test: |r(ixi) - r(2x2) IA(2x2) and 4|r(2x2) - ^(4x4) IA(4x4) along u = 
20, 40, 60, 80. This shows the second order convergence. 

or 

r,v = (54) 
g^u = Right hand side of Equation (55) 

and so on. In any case, ah equations are first order differential equations. Therefore, for example, 
if we apply Equations (|5Q|) and ([5T]) . we can successively solve 

f{Au + u^, Av + Vi) = f{u^, Av + v^) + (RHS of Equation (|32]) at (ii„ Av + Vi)) x Au, (56) 
r(A7x + ii^, + i;^) = r{u^, Av ^ Vi) ^ f{ui, Av ^ Vi) x Au. (57) 

Or, if we use Equations (j52j) and (|53]) . then we can successively solve 

^(A^i + ii„ Av + Vi) = ^(Aii + Ui, Vi) + (RHS of Equation (|33l) at (Aix + u^, Vi}) x Av, (58) 
r{Au -\- Ui^ Av -\- Vi) = r(A'u + ix^, v^) + ^(A'u + ii^, v^) X Av. (59) 

If we use Equations and (j55j), then we can successively solve 

g{Au + Av + v^) = g{u^, Av + v,) + (RHS of Equation 1^ at (ix^, Av + v,)) x Au, (60) 
r{Au -\- Ui, Av -\- Vi) = r(Aix + ii^, v^) + ^(Aix + li^, v^) X Av. (61) 



28 



Although the actual calculations are more difficult than this since many functions are involved, the 
principle is the same. Note that these are based on the first order method (the Euler method). We 
can improve the integration scheme to the second order [26 : to solve the equation 



|=F(.,,), (62) 



we successively solve 



ki = AxxF{xn,yn). (63) 

/C2 = AxxF(xn + Ax/2,7/n + /ci/2), (64) 

= yn^k2^0{{Axf), (65) 

where n is the step and Ax is the step size. It is not difficult to implement the second order 
Runge-Kutta method to the series of equations of this paper. 



Appendix B: Consistency and convergence tests 

In this appendix, we report on the convergence and consistency tests for our simulations. As a 
demonstration, we check the case of hyperbolic symmetry: Vi with Vf = 0.0001, tq = —mo = 100, 
Au — Av — 20, iisheii = '^sheii = 30, St = 0, and = 0.1. 
For consistency, we test one of the constraint functions: 



\9,v - 


-2gd- 






{\9,v\ + 


2gd- 


47rrT„t 





(66) 



around u = 20,40,60,80. Figure [TSl shows that it is less than 1 % except some points, where the 
denominator oscillatory vanishes {g^y ~ 0); this will not be accumulated as one integrates along v. 
Therefore, this shows good consistency. 

For convergence, we compared finer simulations: 1x1,2x2, and 4x4 times finer for around 
u = 20,40,60,80. In Figure [T9l we see that the difference between the 1x1 and 2x2 times finer 
cases is 4 times the difference between the 2x2 and 4x4 times finer cases, and thus our simulation 
converges to second order. The numerical error is < 10~^%. 
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